Static and dynamical dipolar / strain fluctuations in perovskite ferroelectric relaxors 
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We develop a theory to study the characteristics of dipolar / strain fluctuations in perovskite 
relaxors. In addition to the soft TO and TA phonons, we take into account another freedom of 
motion associated with the random hopping of Pb^^ ions between the off-center sites around the 
high symmetry corner site, thus constructing a coupled TO-TA-pseudospin model to describe the 
perovskite relaxors. It is shown that there is a possibility that prior to the on-set of instability of the 
uniform TO-mode (soft-mode), instability of TA mode with nanometer scale modulation (q ~ 0.1a*) 
takes place, which produces static heterogeneous structure concerning polarization as well as shear 
strain. This seems to suggest the intrinsic origin of the heterogeneity in relaxors visualized as 
random distribution of PNR. The phonon spectral density distribution has also been investigated. 
It is shown that when the relaxation time of random hopping of Pb^^ ions is comparable to the TO 
phonon frequency at q ~ 0.1a*, the calculated phonon spectral density reproduces the characteristic 
features called 'waterfall'. 



I. INTRODUCTION 

Recently, the unique properties of ferroelectric relax- 
ors have become one of the central topics in solid state 
physics. The diffuseness of the dielectric response against 
temperature variation has been considered to be due to 
the randomness of the system created by the random oc- 
cupation of B-site ions with different valences. In fact, 
the basic physics of relaxors has been mainly discussed 
in terms of 'random field' to stabilize a glassy state anal- 
ogous to spin glasses. 

One of the intriguing problems to be answered is the 
role of Pb^+ -ions on A-site. In spite of the efforts to 
eliminate Pb ion based on technological reasons, hith- 
erto known perovskite relaxors seem to be restricted to 
Pb compounds, which suggests that at least part of the 
unique properties would be related to Pb ions at A-site. 

From experimental view point, there exist a few unique 
features exhibited by relaxors both in structural as well as 
lattice dynamical aspects providing some 'key concepts' 
to describe characteristic features of relaxors. The key 
concept to discuss the static structure of relaxors is the 
so-called PNR (polar nanoregion). That is, the averaged 
cubic symmetry is locally broken in the temperature re- 
gion where dielectric constant exhibits broad maximum. 
The overall structure is thus visualized as a random dis- 
tribution of PNR embedded on cubic parent phase. 

We notice that there are a few materials which show 
similar intrinsic heterogeneity in the vicinity of phase 
transition point. A group of bcc-based alloys called 
shape memory alloys, which undergo martensitic trans- 
formation, exhibit heterogeneous structure where the 
'embryo's or microdomains of martensite are embedded 
on the austenite (bcc) matrix over a wide temperature 
range. More recently, a group of perovskite manganites 
called CMR substances, which undergo metal-insulator 
transitiou. have been noticed to develop heterogeneous 
structured where microdomains of metallic phase are co- 
existing with insulator phase. The close relationship be- 



tween CMR substances and relaxors was already pointed 
out by Kimura et al.H 

Onukii investigated the origin of the stability of two- 
phase coexistence in alloys and pointed out that the cou- 
pling between the order parameter to the local strain is 
essential to stabilize the heterogeneous structure. Later, 
Yamada and TakakuraQ also arrived at the same conclu- 
sion in the case of CMR substances. 

In this connection, the recent neutron scattering study 
on PMN by Hirota et alJj seems to be very suggestive. 
They carried out the dynamical structure analysis of the 
diffuse scattering and concluded that the displacement 
pattern of each ion in the unit cell contains considerable 
amount of CM (center of mass) non-conserving compo- 
nents. In the language of phonon modes, it means that 
the normal coordinate of the condensing mode is given by 
a linear combination of TO-modc and TA-mode suggest- 
ing the possibility of strong polarization-strain coupling. 

On the other hand, the key concept to characterize 
the dynamical aspect of relaxors is, so-called 'waterfall' 
in phonon spectrum. Gehring et aLH carried out the pio- 
neering neutron scattering study on PZN-8PT and found 
out that the observed high intensity ridge of the neutron 
spectrum did not follow the expected TO phonon dis- 
persion at g < 0.2a*. Instead, it falls down vertically to 
precipitate onto TA dispersion. Similar features are suc- 
cessivefet-pbserved in various relaxors including PMN and 
PZN,Qa'EI indicating that waterfall is indeed the unique 
lattice dynamical characteristic of relaxors. 

Gehring et al0 analyzed the neutron spectrum of PZN 
based on the 'mode coupling' treatment which was uti- 
lized by Harada et ahllil to explain the phonon spectrum 
of BaTiOs. They concluded that in order to reproduce 
the observed spectrum, the TO phonon width should 
show an abrupt change at g ^ 0.1a* where waterfall takes 
place. They claim that the abrupt change is caused by 
the random distribution of PNR in the medium, thereby 
the propagation of the lattice wave with wave length 
longer than the average size of PNR is impeded. While 
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this viewpoint is very attractive, whether the static het- 
erogeneity such as PNR wih cause waterfall type anomaly 
or not is unclear since a static entity will only give rise 
to large momentum transfer of the phonon without caus- 
ing any energy transfer from phonons to other freedom 
of motion and eventually to heat bath. In order to estab- 
lish effective channels of energy transfer, some dynamical 
entity to which phonons are coupled would be needed. 

We consider that the configurational freedom of motios. 
of Pb^+-ion will provide such possibility. It is knowntHi 
that in relaxors the instantaneous equilibrium position 
of Pb^+ ion is slightly shifted from the high symme- 
try corner site and is making random hopping motion 
between the equivalent off-center sites to recover cubic 
symmetry on average. Formally, such freedom of mo- 
tion can be expressed by a stochastic pseudospin vari- 
able. Based on these considerations we construct a suit- 
able model of the relaxor which is characterized by 'cou- 
pled TO-TA-pseudospin system'. In the next section, 
we discuss the TO-TA coupling within the framework 
of 'quasi-harmonic' approximatioii-ivhich was utilized by 
Axe, Harada and Shirane (AllS)tj in order to analyze 
the anomalous TA dispersion in KTaOa. When applied 
to PMN, this treatment suggests the possibility to sta- 
bilize a heterogeneous static structure. In section 3, we 
discuss the neutron scattering spectra including water- 
fall anomaly by taking into account of random hopping 
of Pb^+ ions. The last section is devoted to conclusions 
and discussions. 



II. STATIC POLARIZATION / STRAIN 
FLUCTUATIONS - ORIGIN OF POLAR 
NANOREGION - 

According to the analysis of the diffuse scattering in- 
tensities observed in PMN axound various reciprocal lat- 
tice points by Hirota et al.Jj the structure factor of the 
condensing mode includes substantial amount of atomic 
displacements which does not conserve the center of mass 
(CM) of the unit cell. Since the optical modes at g = 
should satisfy the condition of CM-conservation of the co- 
ordinates, the above observation seems to suggest that in 
PMN, there would be strong coupling between the soft 
optical mode and the transverse acoustic mode around 
g = 0, so that the normal coordinates include consider- 
able amount of CM-nonconserving displacements. In this 
connection, we notice that there are a few perovskite fer- 
roelectric materials, which exhibit anomalous behavior 
in the TA branch upon softening of the TO branch. Fig. 
1 shows typical behavior in the case of KTaOs.Efl It is 
seen that as the TO mode softens, the dispersion of the 
TA branch exhibits anomalous 'concave' curve in a lim- 
ited g-region of 0.02 < q < 0.2. It is noticeable that the 
dispersion shows no 'softening' in the vicinity of the re- 
ciprocal lattice point within q < 0.02. Such anomaly in 
the TA dispersion also indicates the possibility of strong 
coupling between TO- and TA- modes in perovskite fer- 



roelectric systems. 
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FIG. 1: The TO and TA phonon dispersions of KTaOa at 
300K and 20K given in ref. 14. As the TO branch 'softens' 
upon lowering temperature, the TA branch exhibits anoma- 
lous concave curve in a limited g-region of 0.02 < q < 0.25. 
Notice the dispersion shows no anomaly in the direct vicinity 
of the zone center (q — 0) within q < 0.02. 

In 1970, Axe et ahlil discussed the anomalous behav- 
ior of TA branch in KTaOs in the framework of 'quasi- 
harmonic coupling' treatment. They express the dynam- 
ical matrix D(q) of the harmonic potential by taking the 
normal coordinates at g = as the basis functions. The 
off-diagonal matrix elements Dij (q) may then be consid- 
ered as the mode-mode coupling energy between the i'th 
and j'th modes. By taking only TA-TO coupling into 
account the phonon properties (characteristic frequency 
and the corresponding normal coordinates) are explicitly 
given by solving the following secular equation: 



= 0, 



1 : TO 

2 : TA 



(1) 



where ujo(T) is the soft mode frequency at g = which 
is the only temperature dependent quantity, i^y (q) are 
to be expanded in terms of | q | as: 



F,,{q)^f^J\n)q^ + f^f^{K}q\ 



(2) 
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K = q/|q|, (3) 

along a specific direction k. The diagonal elements 
f^j\K)^s may be obtained by comparing with the ex- 
perimental dispersion curves at high temperatures where 

(2) 

no anomaly is observed. In particular, (i^) is directly 
given by the elastic constant for shear strain along k- 
direction. The only essential parameters to be fixed are 
fij^^s with i ^ j describine the coupling between TO- 
and TA- modes. Axe et al.E3 analyzed the experimental 
results on KTaOa and successfully reproduced the TA, as 
well as TO, dispersions at various temperatures by fitting 
these two parameters. 

We further notice that the off-diagonal matrix elements 
of the transformation matrix of the basis functions, S(q), 
become comparable to the diagonal ones in the q-region 
where the TA branch exhibits anomalous behavior. For 
instance, S(q) at q = (0, 1, 0, 0) in KTaOa is given by, 

a^r, 1 n ^^ f ^.880 0.475 \ 

8(0.1,0,0) =(^^0.475 0.880 j- 

This means that in the g-region where anomaly in TA 
dispersion takes place the normal coordinate of TO-mode 
contains large amount of uniform translational displace- 
ments, and vice versa. (TA-mode contains large amount 
of pure optical displacements.) This behavior reminds us 
the characteristics of the dynamical structure factors as 
pointed out by Hirota et al. in the case of PMN. 

It is noticeable that PbTiOa, the-prototype material of 
perovskite relaxors, was reportedll3 to exhibit the same 
type anomaly in TA dispersion. Therefore, it would not 
be unreasonable to assume that perovskite relaxors also 
belong to the materials which experience strong TO-TA 
coupling. Unfortunately, such anomaly can not be proved 
directly by the observation of neutron spectra in relaxors 
because of the extraordinarily large damping character- 
ized by the 'waterfall' phenomena. 

At this stage, we try to construct the 'hypothetical' 
dispersion curves of PMN within the 'quasi-harmonic' 
coupling formalism by completely neglecting the anhar- 
monicity of the potential to cause damping. Among 

(2) 

the parameters defined in eq. (2), (**) is unambigu- 
ously given by the observed elastic constants C44 aad 
l/2(cii-ci2) for /« // [100] and K // [110] respectively.ta 

The remaining fj^j^KYs are assumed to be isotropic (in- 
dependent of k) and to take the same values as those for 
KTaOa given by Axe et al.. 

The results of calculation are given in Fig. 2. The 
dispersion of TA [100] shows the expected concave curve. 
In contrast, the dispersion of TA [110] exhibits a 'dip' 
around q ^ 0.1. Moreover, by subtle adjustment of the 
parameters, the dip becomes more pronounced so that 
wta ^ at g = 0.12. That is, the TA mode has become 
'condensed' to form a modulated static structure with 
modulation period Aq = 27r/go = oTI' 
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FIG. 2: The calculated 'hypothetical' phonon dispersions in 
PMN. The dashed lines are the asymptotic behajvior at g ~ 
determined by the observed elastic constantstj along [100] 
and [110]-directions. The dotted curve corresponds to the 
case when the coupling parameters are slightly modified from 
the values determined for KTaOs. 

The transformation matrix at go — 0.12 is calculated 
to give 

5,0.1,0.12,0)= aS^)- (^) 

Therefore the local structure of the condensing TA mode 
is expressed as the linear combination of uniform trans- 
lational and pure optical displacements with the ratio of 
S21/S22 = 0.50. 

Physically, this implies that the TA [110] mode with 
the wave length of nanometer order, Aq, has the ten- 
dency to become unstable prior to the condensation of 
the uniform (g = 0) TO-mode due to the quasiharmonic 
coupling. Hence the system eventually stabilizes a po- 
larization / strain modulated structure with wave length 
Aq. As has been pointed out, such characteristics of the 
phonon dispersions are not able to be observed by inelas- 
tic neutron scattering experiment obscured by the large 
damping. On the other hand, X- ray diffuse scattering 
intensity, which is proportional to the instantaneous cor- 
relation of fiuctuations, seems to give some indirect in- 
formation on the characteristic frequency w(q). 

The X-ray intensity due to the excitation of the 
phonons belonging to A'th |bxanch is given, irrespective 
of the property of damping,llZI as 

^'(K) = I ^a(K) |^ (6) 

where F^(K) is the dynamical structure factor of the A'th 
mode. In the present system, where TA-mode, rather 
than TO, is assumed to become extremely soft, we may 
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express the X-ray diffuse intensity around each Bragg 
point, K^, in the form: 

^(K) = I FTAiK) I, (7) 

Fta(K) - S22F%{Kh) + S2iF^o{^h), (8) 

where F^^(K/i) and F^Q(Kh) are the dynamical struc- 
ture factors for the pure translational, and pure optical 
displacements at g = respectively. Fig. 3(a) shows the 
calculated X-ray diffuse scattering intensity distributions 
around three Bragg positions in comparison with-the ob- 
served X-ray results by You et al.. (Fig. 3(b))£a While 
overall characteristics of anisotropic distribution are re- 
produced qualitatively, a remarkable discrepancy is seen 
in the diffuse pattern around (400). The observed dis- 
tribution does not show any clear existence of 'satellites' 
which are shown in the calculated contour. This point 
will be discussed later. 

Although it is not shown explicitly in the figure, the rel- 
ative intensities of the diffuse pattern reflects the degree 
of mixing |«f the modes, through the structure factors. Hi- 
rota et al.tl gave the relative displacement for each ions in 
the unit cell which are divided into CM- conserving (TO- 
like) and CM-nonconserving (TA-like) components (See 
eq. (6) in ref. 5). Using the table, the experimental value 
of S21/S22 is obtained as 0.63, which should be compared 
with the calculated value of 0.50 (See eq. (5)). 

III. DYNAMICAL POLARIZATION / STRAIN 
FLUCTUATIONS -ORIGIN OF WATERFALL- 

The neutron spectrum of PMN along [100] in the tem- 
perature region where 'waterfall' phenomenon takes place 
has been analyzed by Gehling et al.Ej based on the cou- 
pled mode treatment utilized by Harada et al. to discuss 
the asymmetric spectra in BaTiOa. In this treatment, 
there are five adjustable parameters to be fitted to the 
observed neutron intensity profile for each q- value. It 
has been shown that, in order to reproduce the entire in- 
tensity distribution in the observed g-range, the effective 
width of TO-modc, Ti{q), changes abruptly at qo where 
'waterfall' takes place. They claim that the abrupt in- 
crease of Fi would be caused by inhomogeneity of the 
lattice due to random distribution of PNR with charac- 
teristic size of '--^ %^ ■ 

While this interpretation is very attractive, the basis 
of the 'coupled mode' formalism on which the treatment 
is based, is a 'homogeneous anharmonic lattice'. The 
damping of phonons are caused by the energy flow to 
the heat bath through anharmonic potential, whence the 
static heterogeneity of the medium would be out of the 
framework of the treatment. 

We take somewhat different standpoint. Besides 
the phonon system, we introduce a stochastic vari- 
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FIG. 3: (a) Calculated X-ray diffuse intensity distributions 
around various types of Bragg reflections, which are compared 
with the observed intensity distributions by You et al.t3 re- 
produced in (b). 



able whose dynamical behavior is only statistically de- 
termined. When such variable is strongly coupled to 
phonons, the life time of phonons would be mainly de- 
termined by the energy flow to the heat bath through 
the random variable. In perovskite relaxors, the config- 
urational freedom of Pb^"*" ions seems to play the role 
of such stochastic variable since it is well establishedtS 
that the instantaneous equilibrium position of Pb^+ ion 
is slightly displaced from the corner site due to covalency 
effect, and it is randomly hopping across the potential 
barrier between the equivalent off-centered sites allowed 
by the averaged cubic symmetry m3m. We may describe 
this freedom of motion of Pb^+ ion by a pseudospin which 
takes on a few distinct values. Thus, our standpoint is 
schematically envisaged as illustrated in Fig. 4. 

As the suitable framework to treat the coupled TO-TA- 
pseudospin system, we utilize Langevin equation which 
describes the motion of the variable under random force 
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FIG. 4: Schematic description of the couphng scheme in TO- 
TA-pseudospin model. Major energy flow from the phonon 
system to heat bath is caused through the coupling to the 
pseudospin system. 



A{t) = nA- J ^{t-t')- A{t')dt + i{t), (9) 



(10) 



(11) 



where A is the state vector whose components are given 
by the independent variables of the system, x is the static 
susceptibihty tensor defined by 



(12) 



In the simple case of the random force with white spec- 
trum, the time development of the averaged value of A(i) 
is given bytj 



with 



(A(t)) = (7-/3)(A), 
7 = [A,A] + (f -f), 



(13) 



(14) 



(15) 



Once the linear equation of motion is established as above 
it is not difficult to obtain the spectral representation of 
the correlation function, (pij{i^) — J {Ai{t)Aj{0))e'^'^'^dt, 
in the formcJ: 



^.,(a;) = [C + ^c./3]^.l + [C-^^/3]^.^ 



(16) 



(17) 



To apply the above general discussions to the present 
system, we define a five-component state vector: 

A+ = (Pi(q),Qi(q),P2(q),Q2(q),a(q)), (18) 

where Pi,Qi are the momentum and the amplitude of 
the TO (i = 1) and TA [i — 2) phonons, and CT(q) is the 
Fourier transformed pseudospin variable: 



a(q) = l/ViV j a^e^^^dr. 



(19) 



As stated above, we assume that the random force, which 
is acting only on the pseudospin variable has the white 
spectrum: 



if fit)) = jS{t). 



(20) 



Physically, 7 corresponds to the relaxation constant of 
the hopping motion of Pb^+-ion. 

The energy of the coupled TO-TA-spin system is ex- 
pressed by, 



H = Y.q{PiicD+u;U<DQii<D) 

K 

-i(P|(q)+c.i(q)Qi(q)) + ij|a(q)| 
+/i2Qi(q)g2(-q) + .gi(3i(q)f7(-q) 
+32(32(q)cr(-q)}, 



.1 



(21) 



where the last three terms give the linear couplings be- 
tween the spin and phonons. 

Using eq. (16), we obtain the phonon spectral density 
in the 5-component system as follows: 



5(q,c.) = 5^ [C + ^^m' + [C- ^^f^, 



(22) 



i=2A 



C = /3-7-/3, (23) 
where (3 and 7 are explicitly given (See Appendix) by. 
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Notice if we eliminate the fifth row and column from /3 
and 7 tensors, the spectral density simply reproduces the 
'hypothetical' phonon dispersions given in Fig. 2 : 



S{uj ± w(q)) 



^^2(q) 



5(^±ci2(q)),(26) 



where w(q)'s are the renormalized TO and TA phonon 
frequencies. In this context, the present treatment may 
be viewed as a natural extension of AHS formalism to 
include the pseudospin freedom of motion. 

In order to visualize general characteristics produced 
by the coupling to pseudospin (stochastic) variable, we 
consider a simpler case of 'single TO phonon-pseudospin 
coupled system.' In this case, we can obtain rather simple 
analytic expression of S{q, uj) as follows 



2„2 



W2(c^2 _ ^2)2 + ^2{^2^2 ^ j,(^2 _ ^2)|: 



r(27) 



(28) 



J' = knT-J 



Similar formula was already given by Yamada et al. O'EHl 
It is worthwhile to notice that the profile of the TO 
phonon spectrum changes drastically as 'two peak — > 
broad single peak — > triple peak' when the relative magni- 
tude of Wo and 7 is changed as 7 ljo ~^ 7 ^ t^o ~* 7 *C 
Wo- (See Figs. 1~3 in ref. 21). Hence, if the relaxation 
constant 7 satisfies the condition: 7 ~ WTo(qo)j the pro- 
file of spectral distribution should exhibit the 'waterfall'- 
like behavior around q = qo • 

Keeping these considerations in mind we use eq. (22) 
to calculate numerically the phonon spectral density dis- 
tribution of PMN along [100]-direction in the region 
< 9 < 0.2. The important parameter 7 has been taken 
as h'y =7meV. The result is depicted in Fig. 5 (a), in 
comparison with the experimental results by Gehring et 
al.. In spite of that we have used only limited number 
of disposable parameters (/12, <7i(= 92), the character- 
istic features of the observed neutron spectrum through 
the whole q-region has been well reproduced. 



(q 0.1a*) takes place, which produces static heteroge- 
neous structure concerning polarization as well as shear 
strain. In the case of PMN, the most probable direction 
of modulation is [110]-direction. This seems to suggest 
the intrinsic origin of the heterogeneity in relaxors visu- 
alized as random distribution of PNR. 

The phonon spectral density distribution has also 
been investigated. It is shown that when the relaxation 
time of random hopping of Pb^^ ions between the 
equivalent off-center sites is comparable to the TO 
phonon frequency at q = 0.1a*, the calculated phonon 
spectral density reproduces the characteristic features 
called 'waterfall' which was observed by neutron spectra 
in PMN and PZN. 

In this paper, we have discussed the possible ori- 
gin of intrinsic heterogeneity in relaxors in terms of TA 
mode instability. Since the energy is simply assumed in a 
quadratic form, the resultant heterogeneity is expressed 
by a harmonic modulation. When higher order terms 
are taken into account, the modulation would become 
more or less 'kink'-like producing well defined boundary 
between the polar and non-polar regions. Recently, Ya- 
mada and Takakura investigated the origin of IC phase 
and two-phase coexistence in Ao.sBo.sMnOs (OMR sub- 
stance) observed around metal-insulator phase boundary 
based on TDGL formalism. Without the higher order 
coupling between 3d-orbital of Mn'^+ ion and local strain, 
the system stabilizes a regular IC structure. However 
when the higher order electron-phonon coupling is taken 
into account, a two-phase coexistent state is stabilized 
in which metallic nanoregion is embedded randomly on 
the parent insulator phase. 

Similar features would be expected in relaxors. Partic- 
ularly, the random distribution of B-site ions provides the 
random pinning center to fix the 'kink' (domain bound- 
ary) position, which will enhance the randomness of the 
spatial heterogeneity. This situation could explain the 
discrepancy between the observed and the calculated dif- 
fuse patterns around (hOO)-type Bragg position as shown 
in Fig. 3. 
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IV. CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have developed a theory to study 
the characteristics of dipolar / strain fluctuations in per- 
ovskite relaxors. In addition to the soft TO and TA 
phonons, we take into account another freedom of motion 
associated with the random hopping of Pb^^ ions be- 
tween the off-center sites around the high symmetry cor- 
ner site, thus constructing a coupled TO-TA-pseudospin 
model to describe the perovskite relaxors. 

It is shown that there is a possibility that prior to the 
on-set of instability of the uniform TO-mode (soft-mode), 
instability of TA mode with nanometer scale modulation 
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APPENDIX A 

We start with general thermodynamical equation of 
motion of a mult-component variable x(= {xi}): 

X = -7X, (Al) 



PZN 




(2,k,0) 



FIG. 5: (a) Calculated phonon spectral density distribuition, 5'(q, tj), along [100]-direction in PMN, which is compared with the 
observed neutron scattering spectrum by Gehring et al.u reproduced in (b). The dashed lines corresponds to the 'hypothetical' 
phonon dispersions when the coupling to the pseudospin is neglected. 



where X is the 'driving force' of the system defined in Using these linearized equations, the spectral representa- 
terms of the thermodynamical potential i^(x) given by tion of the two-time correlation function: 



X 



1 dF 



(A2) 



When F is given in the quadratic form with respect to x 
as, 



a 

X becomes the linear function of x: 



(A3) 



ipa{uj)^ J {x,{t)xi{t + T))e'^'dT (A7) 
is expressed in terms of the coefficients /3, 7 as follows: 



(A8) 



We further notice that when (puiuj) is integrated over uj, 
we obtain the instantaneous (same time) correlation as, 



X = /3 • X. 

Substitution of (A4) into (Al) gives 

X = -(7 • /3)x. 

Conversely, X is expressed as 

X = -C-x 

= (/3-7-/3)x. 



(A4) 



(A5) 



(A6) 



ipii {uj)dLo 



x,{t)xi{t + T))e" 
XiXi) ^ ksTxii 



^drdoj 



(A9) 



where x is the static susceptibility tensor. It is shown 
that within the linear approximation the instantaneous 
correlation is given by. 



(xiXi) = [f3' 



(AlO) 
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In order to apply the above treatment to the present 
case, we consider that the long wave TA and TO phonons 
under consideration may be treated as thermodynami- 
cal variables. Using the energy given in eq. (21), the 
associated thermodynamical potential is expressed in a 
quadratic form; 

+ I a(q) |2 -kT I a(q) ^ +/i2Qi(q)Q2(-q) 



+5iC7(q)Qi(-q) +52CT(q)g2(-q). (All) 



As for the entropy contribution, we have only taken 
into account the configurational randomness of the pseu- 
dospin variables. 

Finally we compare the fundamental equations (13), 
(14), and (15) in the text with (A5), (A9) and (AlO), 
from which the expressions for /3-and 7-tensors are de- 
duced as given in eqs. (24) and (25). 
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